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1. Introduction 



In 1963 Kac et al. [KUH] introduced a statistical mechanics model of particles interacting via 
long, but finite range interactions, i.e. through potentials of the form J 7 (r) = ^f d J{^r), there J is 
some function of bounded support or rapid decrease (the original example was J(r) = e~ r ) and 7 
is a small parameter. These models were introduced as microscopic models for the van der Waals 
theory of the liquid-gas transition. In fact, in the context of these models it proved possible to derive 
in a mathematically rigorous way the van der Waals theory including the Maxwell construction in 
the limit 7 j 0. In mathematical terms, this is stated as the Lebowitz- Penrose theorem [LP]: The 
distribution of the density satisfies in the infinite volume limit a large deviation principle with a 
rate function that, in the limit as 7 tends to zero, converges to the convex hull of the van der Waals 
free energy. For a review of these results, see e.g. the textbook by Thompson [T]. 

Only rather recently there has been a more intense interest in the study of Kac models that 
went beyond the study of the global thermodynamic potentials in the Lebowitz-Penrose limit, 
but that also considers the distribution of local mesoscopic observables. This program has been 
carried out very nicely in the case of the Kac-Ising model in one spatial dimension by Cassandro, 
Orlandi, and Presutti [COP]. A closely related analysis had been performed earlier by Bolthausen 
and Schmock [BS]. These analysis can be seen as a rigorous derivation of a Ginzburg-Landau type 
field theory for these models. Very recently, such an analysis was also carried out in a disordered 
version of the Kac Ising model, the so-called Kac-Hopfield model by Bovier, Gayrard, and Picco 
[BGP1,BGP2]. 

An extension of this work to higher dimensional situations would of course be greatly desirable. 
This turns out to be not trivial and, surprisingly, even very elementary questions about the Kac 
model in d > 2 are unsolved. One of them is the natural conjecture that the critical inverse tem- 
perature (3 c (l) i n the Kac model should converge, as 7 j 0, to the mean-field critical temperature. 
This conjecture can be found e.g. in a recent paper by Cassandro, Marra, and Presutti [CMP]. In 
that paper a lower bound /3 c (t) > 1 + &7 2 | hi7| is proven. A corresponding upper bound is only 
known in a very particular case where reflection positivity can be used [BFS]. 

In addressing this question one soon finds the reason for this unfortunate state of affairs. All 
the powerful modern methods for analyzing the low-temperature phases of statistical mechanics 
models, like low-temperature expansions and the Pirogov-Sinai theory, have been devised in view 
of models with short range (often nearest neighbor) interactions, with possible longer range parts 
treated as some nuisance that can be shown to be quite irrelevant. To deal with the genuinely 
long-range interaction in Kac models, that is to exploit their long range nature, these methods 
require substantial adaptation. The purpose of the present paper is to help to develop adequate 
techniques to deal with this problem - that beyond proving the conjecture of [CMP] will, hopefully, 



also provide a basis for the analysis of disordered Kac models. (Together with possible other means 
not touched by the presented article : most notably with suitably developed expansion techniques 
for long range models). 

The model we consider is defined as follows. We consider a measure space (<S, J-) where 
S = {—1, 1 } z ' ! is equipped with the product topology of the discrete topology on {—1, 1} and T is 
the corresponding finitely generated sigma-algebra. We denote an element of S by a and call it a 
spin-configuration. If A C ZZ d , we denote by a A the restriction of a to A. For any finite volume A 
we define the energy of the configuration a a (given the external configuration a A c ) as 

H 1;A (a A ,a A -) = -\ ^2 J 7 (i - 3)°~i°~j ~ ~ ( L1 ) 

where J 1 {i) = ^y d J(ji) and J : LR d — > IR is a function that satisfies J md dxJ(x) = 1. For simplicity 
we will assume that J has bounded support, but the extension of our proof to more moderate 
assumptions on the decay properties of J is apparently not too difficult. To be completely specific 
we will even choose J(r) = Cdl| x |<i where c<j normalizes the integral of J to one 1 . Here | • | is most 
conveniently chosen as the sup-norm on IR d . 

Finite volume Gibbs measures ("local specifications") are defined as usual as 

< M N = ^_ e -0*,.A(*A,,A.) (1.2) 

J,f3,A 

where ^ A is the usual partition function. Note that under our assumptions on J the local 
specifications for given A depend only on finitely many coordinates of ij. An infinite volume Gibbs 
state jU 7!i g is a probability measure on (S, J 7 ) that satisfies the DLR-equations 

Our first result will be the following 

Theorem 1: Let d > 2. Then there exists a function f(-y) with lim 7 j /(7) = such that 
for all (3 > 1 + /(7), there exist at least two disjoint extremal infinite volume Gibbs states with 
local specifications given by (1.2). Moreover, for 7 small enough, 7(7) < 7(2<j+2)(i+i/d) f or arbitrary 
e > 

Remark: This theorem shows that the conjecture of [CMP] is correct. Together with Theorem 1 
of [CMP] it implies that lim 7 jo/3 c (7) = 1 in the Kac model. While completing this work we have 
received a paper by M. Cassandro and E. Presutti [CP] in which the conjecture of [CMP] is also 

1 The generic name Cd will be used in the sequel for various finite, positive constants that only depend on 
dimension. 



proven, but no explicit estimate on the asymptotics of the function f(-y) is given. Their proof is 
rather different from ours. Although at the moment we make use of the spin flip symmetry of the 
model, the contour language we introduce is also intended as a preparatory step for future use of 
the Pirogov-Sinai theory for non-symmetric long range models. 

We will in fact get more precise information on the infinite volume Gibbs measures in the 
course of the proof. This will be expressed in terms of the distribution of "local magnetization" , 
m x (a), defined on some suitable length scale 1 <C t <C 7 _1 . Given such scale £, we will partition 
the lattice Z d into blocks, denoted by x of side length t. We set Identifying the block x with its 
label x £ we could thus set 

x = {ie% d \ \i-lx\ <£/2} (1.4) 

We then define 

In the sequel we will assume that all finite volumes we consider are compatible with these blocks, 
that is are decomposable into them. We will also assume that j£ is an integer. For any volume A 
compatible with the block structure, we denote by Ma C T\ the sigma-algebra generated by the 
family of variables {m x (a)} xe A- The block variables will be instrumental in the proof of Theorem 
1. However, they are also the natural variables to characterize the nature of typical configurations 
w.r.t. the Gibbs measure. We should note that this first step of passing to the variables m x (a) is 
also used in [CP], in fact it is used in virtually all work on the Kac model. 

The remainder of this article is organized as follows. In Section 2 the distribution of the block 
spins are formally introduced and the block-spin approximation of the Hamiltonian is discussed. 
In Section 3 we introduce our notion of Peierls contours and prove our theorem through variant of 
the Peierls argument [P]. 

Aknowledgements: We thank Errico Presutti and Marcio Cassandro for sending us a copy of 
their paper [CP] prior to publication. M. Zahradm'k also acknowledges useful discussions with E. 
Presutti on Kac models in general and about their recent preprint in particular. We would like to 
thank also the home institutions of the authors and the Erwin Schrddinger Institute in Vienna for 
hospitality that made this collaboration possible. 



2. Block spin approximation 

All the questions we want to answer in our model will after all concern the probabilities of 
events that are elements of the sigma- algebras My for finite volumes V. If A € A4y is such an 
event and A D V, we have the following useful identity 

</3,a(-4) = Y ^7,/3,a( ct a\v)^,v' 7AC (^) 

CT A\V 

cr A \ v m x ,xeV 

The sum over m x runs of course over the values {— 1, — 1 + 2£~ d , . . . , 1 — 2£~ d , 1} Note that we may, 
if J has compact support, assume without loss of generality that A is sufficiently large so that the 
local specification ^ A ^ VAC does not depend on rj. We will therefore drop the r\ in this expression. 

The main point which makes the Kac-model special, is that the Hamiltonian is "close" to a 
function of the block spins. Namely, we may write 

H 1 y{a v ,a V o) = -\ ^ ^ J 1 (i,j)a i a j - ^ ^ J J (i,j)a i a j 

x,y(zV i£x,j£y x£V,y£V c i€x,j£y 

= -1 E M e ( x ~y)) Y Gi<7 i 

x,y£V iex,jey 

- £ M£(x-y)) £ aw (2.2) 
xev,yev c iexjey 

~^Y Y i^-j) - Mtfr-vWwj 

x,yeV i€x,j€y 

= H j},v( m v(^v), , m V o(a V c)) + AH^ t£y (a v ,a V c) 
where we have set (recall that J^ilx) = £~ d Ji 1 (x)) 

H( y },v( m v,in V c) = -£ d \ Y J 7 e(x - y)m x m y - £ d ^ J 1 n(x - y)m x m y (2.3) 

x,yeV x€V,yeV c 

and 



AH lAV (a v ,a V o) = -\ ^ ^ [J^i - j) - J^{£(x - y^a^j 

x,y€V i€x,j€y 

- Y Y [ J -^ ~ ft ~ J ^ x - y ^ Ui °i 



(2.4) 



x€V,yeV c iex,j€y 

Lemma 2.1: For any V C Z d , 

snp\AH lAV (a v ,a V c)\ < c d j£\V\ (2.5) 
where d is some numerical constant that depends only on the dimension d. 



Proof: This fact is well-known and simple for all Kac models. In our case it follows from the 
observation that [J 1 {i — j) — J 1 {l{x — y))\ = 0, unless \x — y\ ~ l/(j£).Q 

As consequence of Lemma 2.1 we get the following useful upper and lower bounds for the 
distribution of the block spins: 

-f3t d H { °\ v {m v ,m v c) yi n? IT 

Of course 



2~ e (i±L ed ), if£ d /m x /2e% (2 _ 7) 
0, else 



and thus, by Sterling's formula, 



r "'^) =e "" /K)+ ° (M) (2 - 8) 



where /(m), for m E [—1, 1] is 



2 -J- 177, \ — jyi 

I(m) = ln(l + m) + ln(l - m) (2.9) 



Therefore we define 



^ 7 ,W,v( m v. m vO = -5 J-ye(x-y)m x m y - ^ J ie (x-y)m x m y +(3 1 ^ 7(m x ) (2.10) 
to get 

Lemma 2.2: For any finite volume V and any configuration my, we have 

< p-/K' i S 1 ,)J,<,v(mv,mvc(ffyc)) 
^j,/3,vK m V) y p -l3^E 1<f3 ^ v {m v ,m vc (a V c)) e ^- L1 > 



Remark: I will be chosen as tending to infinity as 7 tends to zero. The idea is that that E^^^y is 
in a sense a "rate function"; that is to say, E 1 ^^y alone determines the measure since the residual 
entropy is only of the order ^4^|V|. The problem is that this is only meaningful when we consider 
events A for which the minimal E^^^y is of order |V| above the ground state to make sure that 
neither the residual entropy nor the error terms in (2.11) may invalidate the result. We will have 
to work in the next section to define such events. 

It is instructive to rewrite the functional E lt p^y in a slightly different form using that 
—m x m y = \(m x — m y ) 2 — \{m x + m y ) (we drop the indices 7,/?, I henceforth but keep this 



dependence in mind). We set 



E v (m v ,m V c) = \ ^ J lt (x-y){ 




m x - m y ) 



x,y(HV 



x£V,yeV c 



(2.12) 



x£V 



where fp is the well-known free energy function of the Curie- Weiss model, 




(2.13) 



Then 



Ev{mv,myc) = Ey{my ,myc) — Cy(myc) 



(2.14) 



where 




(2.15) 



xev,yev c 



depends only on the variables on V c . 

The form Ey makes nicely evident the fact that the energy functional favours configurations 
that are constant and close to the minima of the Curie- Weiss function fp(m). 



3. Peierls contours 

In this Section we define an appropriate notion of Peierls-contours in our model and use this to 
proof Theorem 1 by a version of the Peierls argument 2 . The general spirit behind the definition of 
Peierls contours can be loosely characterized as follows: We want to define a family of local events 
that have the property that at least one of them has to occur, if the effect of boundary conditions 
does not propagate to the interior of the system. Then one must show that the probability that 
any of these events occurs is small. We will define such events in terms of the block spin variables 
m x (a). More precisely, since it is crucial for us to exploit that the new interaction is still long 
range 3 , contours will be defined in terms of the local averages, cp x (m), and the local variances, 
tp x (m), defined through 



While the proof of [CP] is also based on a Peierls argument, their definition of Peierls contours is completely 
different from ours. 

3 For that reason it is not possible to directly use the methods developed in [DZ] for studying low temperature 
phases of short range continuous spin models, although some of the ideas in that paper are used in our proof. 




(3-1) 



y 



Then define the sets 



y 



(3.2) 



r={x|||<k(m)|-m*(/3)|>Cm*(/3) or ^(m) > (C™*(/?)) 2 } (3.3) 

where m* (j3) is the largest solution of the equation x = tanh /3x, that is the location of the non- 
negative minimum of the function fp. We recall (see e.g [E]) that m*((3) = if (3 < 1, m*{(3) > if 
(3 > 1, lim^ioo m*((3) = 1 and lim^i ^grr^ = 1- To simplify notation we will write m* = m*(/3) in 
the sequel. C, £ < 1 will be chosen in a suitable way later. Note that if the boundary conditions are 
such that say <f) x {m(rj)) « +m*, then, if the configuration near the origin is such that (j) (m(a) < 0, 
there must be a region enclosing the origin on which <\> takes the value zero and thus belongs to 
£. For a reason that will become clear later, in a first step we will regularize this set. For this we 
introduce a second blocking of the lattice, this time on the scale of the range of the interaction. 
The points u of this lattice are identified with the blocks 

u = \x G Z d \ \x - u/(-/£)\ < l/(2 7 £)} (3.4) 

just as in (1.4). We write in a natural way u(x) for the label of the unique block that contains x. 
We will call sets that are unions of such blocks u regular sets. We put 

r = jx|u(x) nf / 0} (3.5) 

For some positive integer n > 1 to be chosen later, we now set 

T = {x\ dist(x,£ ) < n^fT 1 } (3.6) 

where dist is the metric induced by the sup-norm on IR d . n will depend on (5 and diverge as j3 | 1. 
The precise value of n will be specified later in (3.48). Notice that this definition assures that the 
set r is a regular set in the sense defined above.. Connected components of the set set T together 
with the specification of the values of m x , x £ T are called contours and are denoted by T. For 
such a contour, we introduce the notion of its boundary dT, in the following sense: 

dT = {x G £ I dist(x,£ c ) < n( 7 £)" 1 } (3.7) 

Note that by our definition of T we are assured that dT n £ = 0. We denote by 

D ± = {x I |^(m)Tm*|<Cm*}nr c (3.8) 

and call these regions i-correct. Each connected component of the boundary of T connects either 
to D + or D~ . We will denote such connected components by dTf and dT~ , respectively. 



For a connected set F we denote by int F the simply connected set obtained by "filling up the 
holes" of F. This set is called the interior of a contour. The boundary of int F will be referred to 
as the exterior boundary of F. The connected component of dF that is also the boundary of int F 
will be called exterior boundary of F and denoted by dF ext . 

The strategy to prove Theorem 1 is the usual one. First we observe that if boundary conditions 
are strongly plus, then in order to have that, say, \(j)o{m) — m*\ > (m*, it must be true that there 
exists a contour F such that <G intT . Thus it suffices to prove that the probability of contours 
is sufficiently small. This will require a lower bound on the energy of any configuration compatible 
with the existence of F, and an upper bound on a carefully chosen reference configuration in which 
the contour is absent. We will show later (Lemma 3.8) that a lower bound on the energy can easily 
be given in terms of the functions (f> and ip, a fact that motivates the definition of F. The long 
range nature and of the interaction and the fact that the m x are essentially continuous variables 
require the construction of the extensive "safety belts" around this set in order to assure an effective 
decoupling of the core of a contour from its exterior. The crucial reason for the definition of contours 
through the nonlocal functions 4> and ip is however the fact that these are "slowly varying" functions 
of x for any configuration m. Therefore, even if the core F is very "thin" (e.g. a single point), 
one can show that on a much larger set \(f>o(m) — m* \ or ip x {m) must still be quite large (e.g. half 
of what is asked for in F) . This guarantees that in spite of the very thick "safety belts" we must 
construct around F, the energy of a contour compares nicely with its volume |r|. 

We will now establish the "decoupling" properties. For this we must establish some properties 
of the configuration m on dF that minimizes Eqy for given boundary conditions. 

Definition 3.1: A configuration m^ 1 is called optimal if m opt minimizes Ey(my ,myc) for a 
given configuration myc . 

An important point is that away from F, due to our definition of contours configurations must 
be close to constant in the following sense: 

Lemma 3.2: Assume that dist(x,F) > 1/(7^). Then 




(3.9) 



y 



and 



(ii) for any V C F 




(3.10) 



where the sign depends on whether x (m) is positive or negative in the region. 



Proof: The proof of (3.9) is straightforward from the definition of £ in (3.3) and (3.10) follows 
from (3.9) by the Schwartz inequality. <0 

We will now establish properties of an optimal configuration on regular sets with boundary 
conditions that satisfy properties (3.9) and (3.10). 

Lemma 3.3: Let V be a regular set. Then there exists Q > depending only on the dimension 
d such that if myo is a boundary conditions of + type for which (3.9) and (3.10) hold with ( < (d, 
then for all x € V , \m° pt — m*\ < m*/2. The corresponding statement holds for — type boundary 
conditions. 

Proof: We see from (2.10) that we must have 4 for y € V 

= -^—E v (m v , m V c) = IT 1 1' (my) - <p y {m) (3.11) 

till I iy 

(3.11) can be written as 

m y = tanh ((3<t> y (m)) (3.12) 

We may tacitly assume that 4> y (m) is positive (this assumption will be shown to be consistent). 
Since m* is a stable fixpoint of the function tanh (im that attracts all points on the positive half 
line, it follows that | tanh ((3(/) y (m)) — m*\ < \(j> y (m) — m*\ and in particular, if 4> y (m) < m*, 
tanh ((3cj) y (m)) > 4> y (m), while for 4> y (m) > m*, tanh ((3(f>y(m)) < cf> y (m). We will first show that 
m° pt > m*/2. Let x € V denote a point where 

m x = inf {m y \m y < m*} (3.13) 

yev 

If m x = m*, there is nothing to proof. But if m x < m* , then (3.13) can only be satisfied if 
dist(x,<9V) < 1/(7^). For such points we can write 



m x - m* > ^2 J lt {x - y)(m y - m*) + ^ J lt (x - y)(m y - m*) 

y£V y€V c 



> (m x - m*) ^2 J ie( x ~ v) ~ 2 C m * J ^ J ^ x ~ y ^ 
yev V y ev c 



(3.14) 



where the second line follows by (3.10). Hence 



m x — m > ; = (3.15) 

Y,yev- Jje(x-y) 



4 We ignore the fact that m x takes only discrete values and look for the optimal solution in the space of real- 
valued m. The point is that given such a solution, a discrete valued approximation can be constructed that differs 
in energy by less than |r|/£ d which is negligibly small. 



On the other hand, (3.14) holds for any other point y G V as well, and inserting this into the first 
line of (3.14) we get 



m x -m* > (m x - m*) Y Y J lt (x - y)J 1 e{y - z) - AC,m* (3.16) 

y€V z€V 



Clearly we have won if either 



l-^2^2J ie (x-y)J^(y- z)>8( (3.17) 



or 



Y M x ~y)>^c (3.18) 

V y^ Vc 

Due to the fact that V is composed of cubes of sidelength of the range of the interaction, this 
follows from simple considerations if ( is smaller than some dimension dependent constant. (Here 
is the reason for our definition of £ ). In fact, 

i - Y Y j t £ ( x ~ y) J i^y - z ) = Y J i^ x ~y) + Y Y J ^ x ~ y) J ^(y - z ) (3-i9) 

y£V zeV y€V c y€V zeV c 

The point is that the second term on the right hand side of (3.19) cannot be too small as long as 
dist(x, V c ) < 1/(7^), for regular V (if V is not regular, this statement does not hold, of course; 
just consider a thin long spike entering into V and let x be near the tip of the spike!). In fact, the 
worst situation here occurs if x is at a distance r/("f£) from a "corner" of V c . One easily verifies 
that even in this case 



Y Y J ^ x ~ y)' 1 ^ -*)> 2" (d+1) f ds (r + s) d -\l - s) d 

y ev zev^ ^° 

> 2"( d+1 ) f ds s d -\l - s) d = 2"( d + 2 )- ((d " 1)1)2 
Jo 



(2d - 1)! 

so that (3.18) is verified if 4£ is smaller than this number. The numerical value of that bound can 
of course be improved, but we do not seek to do that. 

Having established that m x > m*/2 in 7, a trivial computation shows that our starting 
assumption that (f) x {m) > is also verified. Thus we have proven that m x pt > m* /2. In the same 
way one shows also that m° pt < 3m* /2 which concludes the proof of the lemma. 

In the sequel the notion of n-layer set defined in the following definition will be convenient. 

Definition 3.4: A regular set V is called a n-layer annulus, if there it is of the form 

V = jx G V c I dist(x, V) < n( 7 ^)" 1 } (3.21) 



for some connected set V that is composed of blocks u. The sets 

V k = {x g V c | (k - l)(7f) -1 < dist(x,V) < fc(7£) -1 } (3.22) 
are called the k — th layers of V. 

Note that the sets dT are by their definition n-layer sets. 

We are interested in some properties of optimal configurations on n-layer sets. For this we will 
use the following simple fact about the function fp, that may be found e.g. in [BG] 

Lemma 3.5: Let fp(m) = f3~ l I(m) - \m 2 . Then, for all m G [-1, 1] 

U(m) - fp{m*) > c(0) (H - m*) 2 (3.23) 

where 

m s Jn^, _ 1 (3 M) 
has the property that c((3) > for all (3 > 1 and 

Pli (m*) 2 12 v ' 

From this we will derive the following Lemma (The analog of this Lemma for short range and 
purely quadratic Hamiltonians appeared already in [DZ]). 

Lemma 3.6: Let V be an n-layer set with n > r/c(j3). Then there exists a layer Vk in V such 
that 

£ (m^) 2 < 2-^(m*) 2 (|y 1 | + \V n \) (3.26) 

x€V k 



Proof: Let us set u x = \m x \ — m* and and use the abbreviation 

x€V k 

and analogously for other functions. Then it is obvious from (2.12) that for any configuration, 

n-l 

E v \ Vl \ V2 (m v \ Vl \ Vn ,m VlUVn ) >J2c((3)\\u Vk \\l+ £ fp(m*) (3.28) 

k=2 xeV\Vi_\V n 

On the other hand, we may consider a configuration that equals m opt on V\ and V n and has 
m x = m* for all x € V\Vi\V^. For this configuration 



E v \ Vl \ V2 (m v \ Vl \ Vn = m*,mgj V J = - £ J l£ (x - y) (m° pt - m*) 2 + £ fp{™*) 

(3.29) 



2 

yGVxUVn 



By the definition of m opt , it must thus be true that 

> £y\Vi\v 2 ( m y\Vi\v„> m ^uv-J - E v \ Vl \ V2 (m°^ Vi 

n-l 



opt * opt 

\v„ ~~ m ' m yiuv„ 



>^c{l3)\\u Vk \\l- - J^e(x - y) (m° y pt - m*)' 



k=2 
n-l 



(3.30) 



>2>O0)|K 



2 

fcll2 



fc=2 



opt 1 1 2 
V-lll2 



+ u 



opt m2\ 



Thus, for any q < n/2, we have 



9+1 
?c(/3) inf L 

from where 



opt 1 1 2 

V fc I12 



II 9 "I - ||lty 



opt 



n+l-k ' 



n-l 

fc=2 



I opt,.2 , 11 opt 1 

\ u V k H2 + ll u v„+i-J 



< 



\u 



opt 1 1 2 1 In opt 1 1 2 



V'i 



; + 



V'i 



9+1 

inf 

A: = 2 



K UI«vJll + ll«v n+1 _J|l] < 



1 



opt 11 2 



2gc(/?) "2 + IK, H2j 

If q is chosen as the smallest integer greater than l/c(/3) this shows that there exist 2 < k < q + 1 
such that 



opt || 2 



(3.31) 



(3.32) 



l«vj| 2 + \\ u v n+1 



M < 



opt 1 1 2 
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Iterating this construction, and using that by Lemma 3.3 



+ ll« : 



opt 1 1 2 

Vi H2 



1 



u 



t\\2+h\\<:\\t<^™*r(\vi\+\v n \) 



1 ||„,opti|2 



1 



(3.33) 
(3.34) 



we arrive at the assertion of the lemma. 



We are now ready to construct our reference configuration and give an upper bound on its 
energy. For given contour T and compatible external configuration m on T c we call m opt the 
configuration on T that minimizes the energy for a given core T. Clearly such a configuration is 
also an optimal configuration on dT in the sense of Definition 3.1. Thus by Lemma 3.2 we know 
that in each connected component dTf of the boundary of T there exists a layer Cf of thickness 
1/(7^) in drf such that \\m°f ^m*\\ 2 2 < 2- r \(m*) 2 \[\Vi{dTf)\ + \V n {dYf)\} For given Cf we 
decompose drf into the two sets 

drf tin = {x e dTf\Cf |,dist(x, J D ± ) > dist(£f ,£>*)} (3.35) 

and 

dT± out ee dVt\dVf. n (3.36) 



Without loss of generality we assume that the exterior boundary of our contour is attached to 
the +-correct region. We now define the reference configuration m re f 



m°J>\ 

-m%*, ifxedAr out 

m* , for all other x € T 

m x , for x € D + 

—m T < for x £ D~ 



if x € dAf t 



m r x ef = < 



(3.37) 



Lemma 3.7: Let m re f be defined in (3.37). Then for any compatible external configuration we 
have that 

E L (rn? f ,mff)<Y,E drtut (m^ , m r ) + £ 2-^(m*) 2 [|y 1 (drf )| + \V n (dVf )|] 

a;Gr\ar out 

(3.38) 



Proof: The proof of this estimate is obvious from the definition of m re f and Lemma 3.6. Note 
that in the terms E gr ± ( m dr ± , mr= ) the interaction energy between dTf out and <9r^ in is not 

i,out \ i,out / 

counted. <0 

Of course the configuration m re f does not contain the contour T. It remains to find a lower 
bound on the energy of any configuration m that does contain a contour with given V. 

To do this, we use the following observation. 

Lemma 3.8: Let U,V,W C % d be any three disjoint sets such that for all y G U U W, 

ExGC/uwuv -M 21 - y) = !■ arid /° r anyyeU J2 xe uuw J i^ x ~ v) = 1 T/ien 

#vuE/uw(WuE/uw,m(vut/uwo c ) > 3 £ ^x(^) + I £ [//?(^x) + //3 (<i>x(m))] 

X€U xeUUW (3.39) 

xev 



Proof: The proof of this lemma is a simple, but, mainly because of boundary effects, somewhat 
lengthy computation that we do not wish to reproduce here. To get the idea, note that in infinite 
volume we have (formally) 

- \ m x m y J ll (x -y) + (3' 1 £ I(m x ) 



E 



!£m x x (m) +P 1 £/(m a .) 

X X 



(3.40) 



+ -fJK) + -/rV*(/(m)) 



where we have put (f) x (I(m)) = Yl y Jje( x ~ y)I( m y)- The last line is obtained by inserting the 
identity 1 = J 1 t(x — y) in the I(m) term and changing the order of summation in the resulting 
double sum. Using the same trick for the first term in the last line of (3.40), and using that, since 
/ is a convex function, (f> x (L(m)) > L(4> x (m), one gets that 



E 



(3.41) 



is a lower bound for (3.40). Trying to repeat this computation in finite volume and carefully dealing 
with the boundary terms leads to the more complicated looking formula (3.39). 

The main point in the estimate (3.39) is that it allows to bound the energy of a configuration 
from below in terms of (f) x {m) and ip x {m) alone. Namely, taking for V and U U W the layers Cf 
and the regions "within" Cf, we see that for any configuration 



E r (m r ,m r c) >y^E Qr ± (m ar ± ,m r c) 

— — — * * i,out \ i,out — / 



i,± 

+ \ £ [U{^{m))-U{m*)] + \ £ Mm) 

xer\dr out *er\er out K 6 - ) 

dist(x,ar „ t )>i/( 7 <?) 

+ £ Mm*) 
xer\ar out 

Next we show that bot 4> x (m) and tp x (m) has nice continuity properties. 

Lemma 3.9: There exists a finite constant Cd depending only on the dimension d such that for 
any contour V, ifT denotes the set 

£={» I (3-43) 

then for allyef, \\<f> y (m)\ -m*\ > or ^(m) > ^p^. 

Proof: Since \m x \ < 1, it is a simple geometric fact that 

\4>x(m) -<f> v (m)\ < c d \x-y\j£ (3.44) 

and 

|^ x (m) - ^ y (m)\ < Ac d \x - y\j£ (3.45) 

for some geometry dependent constant a. Since on T, \4> x (m)\ > (m* or ip x (m) > ((m*) 2 , the 
assertion of the lemma follows. <0 

Remark: The estimates of Lemma 3.9 are very crude. We expect that they can be improved 
considerably. 

A further simple geometric consideration shows on the other hand that T cannot be too small 
compared to T, namely 

Lemma 3.10: There exists a numerical constant c' d depending only on the dimension d such 
that for any contour V , we have that 

El^t^kl (3-46) 



I r l 

Proof: Note that W- is maximal if T consists of a single point in which case (3.46) is obvious. 

Combining the upper bound on the energy of m re f from Lemma 3.7 with the lower bound 
(3.42) obtained from Lemma 3.8 applied for the optimal configuration, using the fact that that E 
and E differ only by a constant that depends only on boundary conditions, and finally employing 
Lemma 3.10 we arrive at 

Proposition 3.11: Let T = (F, m) be a contour with fixed T. Then there exists a reference 
configuration m re ^ in which T does not occur such that 

E L (m L ,m r )-E L (m^,m?j) > \F\ - |(m*) 2 2 -^) |£| (3.47) 

where c d is a finite dimension- dependent constant and c(/3) is the constant from (3.24). 

Proof: We bound £r(%,% c ) from below by the the the corresponding energy of the configura- 
tion m of lowest energy for given T; on the belt of the contour this provides a optimal configuration 
in the sense of Definition 3.1. The same configuration is used in the construction of m re f . After 
the obvious cancelations and using (3.46) and the fact that c(f3) < 1, we get the assertion of the 
proposition. <0 

We must now begin to choose our parameters. We want the Peierls condition, i.e. that the 
coefficient of |£| in (3.47) is positive and as large as possible. The most convenient choice appears 
to choose n in such a way that 

2 _ nc(0) _i cmcm*) 2d 

2 c d (n+l) d 

Calling the solution 5 of this equation n*, we get the Peierls estimate 

W ( \ W ( ref ref\ . 1 4(3) ((m* ) 2d+2 

It is not difficult to verify that 



(3.48) 



n*<C^-ln 

c(/3) 



c(/3)(Cm 



*\2d 



2c d 

for some numerical constant C, if c((3) is sufficiently small. 



(3.50) 



This estimate on the energy difference will only be useful for us if it is large compared to the 
error terms arising from the block approximation. That is to say, we must make sure that 

l«(ffl({m-)»« 



16 c d (n+l) d 



5 By this we will of course understand the smallest integer larger than or equal to the "real" solution 



(the two Cd in this formula are a priori not the same quantities). This gives a relation between 
temperature dependent quantities on the one hand and j£ on the other. It does not impose any 
choice on the parameter £. This arises from the last condition, the comparison between the energy 
of a contour and the entropy, i.e. the number of configurations m on F and of shapes F with fixed 
volume |r|. Even the crudest estimate shows that this number is smaller than £ d ILI(7 d IL1 so that 
(3.51) is complemented by the condition 

1 c(/3)(Cm*) 2d+2 



f3£ a 



- c d -y£ 



>d\nl + \nC (3.52) 



16 c d (n + l) d 

which requires £ to be sufficiently large. In fact we may choose £ as 

£ " 7 c d {n + iy (3 - 53) 

which inserted into (3.52) gives the final condition of (3 in terms of 7. It is clear that for any (3 > 1, 
i.e. c(/3) > and m* > 0, this condition can be verified by choosing 7 sufficiently small. Thus 
using Lemma 2.2 we proved the analog of the Peierl's argument here, namely that the probability 
of a given contour F is smaller than exp(— c|r|| ln£|) which in turn implies that the probability that 
the origin is in the interior of a contour is close to zero (in fact of the order exp(— c(3n d \ In €| )) . 
Moreover, by inserting the asymptotic behaviour of m* and c((3), one verifies easily that if we put 

(3-1 = 7(2d+2Ki+i/d) (3.54) 

for arbitrary e > 0, then (3.52) is verified when 7 is sufficiently small. This gives thus the claimed 
bound on the behaviour of the critical temperature as 7 tends to 0. 

This concludes the proof of Theorem 1.0 "0 

Remark: Let us recall some consequences of what we have just proven: if V denotes the union of 
the interiors of all the contours of a given configuration than the Gibbs probability of the event 

dist(i,V c )>r (3.55) 

is independent of the choice of the point i £ Z d and behaves like exp(— Cr) where C = C(/3,7). 
This implies for example the following statement: The probability of the event that the support 
of all contours surrounding a given point is infinite is equal to zero. One could even refine such 
a statement, giving a more precise meaning to the intuitive idea that "almost all configurations 
(of the mesoscopic observables m) in the translation invariant + Gibbs ensemble have their local 
averages (in the sense of the variables (f> x (m)) in the vicinity of m* except of some (rare, but 
uniformly distributed throughout the lattice) "islands" . (This is the appropriate rephrasing of the 
statement in Sinai's book [S]). 



References 

[BGP1] A. Bovier, V. Gayrard, and P. Picco, "Large deviation principles for the Hopfield model and 
the Kac-Hopfield model", Prob. Theor. Rel. Fields 101, 511-546 (1995). 

[BGP2] A. Bovier, V. Gayrard, and P. Picco, "Distribution of overlap profiles in the one dimensional 
Kac-Hopfield model", WIAS-preprint 221, submitted to Commun. Math. Phys. (1996). 

[BFS] J. Bricmont, J.R. Fontaine, and E. Speer, "Perturbation about the mean field critical point", 
Commun. Math. Phys. 86, 337-362 (1982). 

[BS] E. Bolthausen and U. Schmock, "Convergence of path measures arising from a mean field or 
polaron type interaction", Prob. Theor. Rel. Fields 95, 283-310 (1993). 

[CMP] M. Cassandro, R. Marra, and E. Presutti, "Corrections to the critical temperature in 2d Ising 
systems with Kac potentials", J. Stat. Phys. 78, 1131-1138 (1995). 

[COP] M. Cassandro, E. Orlandi, and E. Presutti, "Interfaces and typical Gibbs configurations for 
one-dimensional Kac potentials", Prob. Theor. Rel. Fields 96, 57-96 (1993). 

[CP] M. Cassandro and E. Presutti, "Phase transitions in Ising systems with long but finite range", 
preprint (1996). 

[DZ] R.L. Dobrushin and M. Zahradmk, "Phase diagrams for continuous-spin models: an exten- 
sion of the Pirogov-Sinai theory", in "Mathematical Problems of Statistical Mechanics", R.L. 
Dobrushin, Ed. (D. Reidel Publishing Company, Dordrecht, 1986). 

[E] R.S. Ellis, "Entropy, large deviations, and statistical mechanics", (Springer- Verlag, Berlin, 
1985). 

[KUH] M. Kac, G. Uhlenbeck, and P.C. Hemmer, "On the van der Waals theory of liquid-vapour 
equilibrium. I. Discussion of a one-dimensional model", J. Math. Phys. 4, 216-228 (1963);'TI. 
Discussion of the distribution functions", J. Math. Phys. 4, 229-247 (1963); "III. Discussion 
of the critical region", J. Math. Phys. 5, 60-74 (1964). 

[LP] J. Lebowitz and O. Penrose, "Rigorous treatment of the van der Waals Maxwell theory of the 
liquid- vapour transition", J. Math. Phys. 7, 98-113 (1966). 

[P] R. Peierls, "On the Ising model of ferromagnetism" , Proc. Cambridge Phil. Soc. 32, 477-481 
(1936). 

[S] Ya.G. Sinai, "Theory of phase transitions. Rigorous results", (Pergamon Press, New York, 
1982) 



[T] C.J. Thompson, "Classical equilibrium statistical mechanics" , (Clarendon Press, Oxford, 1988). | 



